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Abstract 

The differential pt spectrum for vector boson production is computed at next-to- 
leading fixed order and including the resummation of threshold logarithms at next-to- 
next-to-leading logarithmic accuracy. A comparison is made to atlas data on direct 
photon and W production at high transverse momentum p?, finding excellent agreement. 
The resummation is achieved by factorizing contributions associated with different scales 
using Soft-Collinear Effective Theory. Each part is then calculated perturbatively and 
the individual contributions are combined using renormalization group methods. A key 
advantage of the effective theory framework is that it indicates a set of natural scale 
choices, in contrast to the fixed-order calculation. Resummation of logarithms of ratios 
of these scales leads to better agreement with data and reduced theoretical uncertainties. 



1 Introduction 



Cross sections for vector boson production are relatively clean observables at hadron collid- 
ers. The differential px spectra for photons, W bosons and Z bosons in particular, provide 
excellent benchmarks to test the standard model as well as to measure parton-distribution 
functions (PDFs). An important application of such measurements is to compare and validate 
different precision calculations, performed at fixed order or including resummation. Vector 
boson production, therefore, gives us a rare handle to gauge the importance of higher order 
perturbative effects and power corrections. In this paper, we perform such a comparison on 
the high px photon |T] and W boson spectra [2] measured by the atlas collaboration at the 
LHC, using around 35 pb _1 data. 

Direct (or prompt) photon production is the production of a hard photon in association 
with a jet. The cleanest direct photon observable is the inclusive photon px spectrum, which 
can be measured independently of any jet definition. At low px, there is a large background 
from 7T° and other hadronic decays, which are often corrected for by demanding that the photon 
be isolated. In the ATLAS study [1] , the isolation criteria was that there should be no radiation 
with less than 4 GeV of energy in a cone of radius R = 0.4 around the photon. An advantage 
of studying the direct photon spectrum at high pt is that there is little background of a hard 
photon coming from background processes and isolation becomes unnecessary. Formally, the 
backgrounds provide only power corrections in this region. 

W and Z production have smaller cross sections than photons, especially after paying 
the cost of a branching ratio to leptons, but do not require isolation. The W spectrum is 
particularly challenging to measure since it requires an understanding of the missing energy, 
which unlike the lepton px, requires mastery of systematic effects over the entire detector. 

The photon, W and Z production rates have been known at the next-to-leading order 
(NLO) for some time IH El HH 13 IB] • In this paper, we take leading order (LO) to refer to 
the leading order in which the vector boson V has non-zero px- So this is a tree-level 2—7-2 
scattering process. NLO is one order beyond this, which includes 1-loop corrections to the 
2—7-2 processes as well as 2 — y 3 real emission graphs. While the inclusive V production rates 
are known at NNLO, the differential px spectra are only known at NLO. These corrections 
are implemented in Monte-Carlo integration programs to provide the NLO distributions, such 
as qt 0, mcfm [TO] , fewz [TTJ [T2j and DYNNLO [T3"] . 

Beyond NLO, the theoretical calculation of the vector boson spectrum is extremely chal- 
lenging, and the NNLO result is not yet known. In the absence of this result, one can improve 
on NLO by adding in partial results at higher orders. In some cases, such as at low pt, 
this is absolutely critical. The fixed-order calculation diverges at small px so one needs to 
resum logarithms of the form ln(pr/my) to get even qualitative agreement with data. The 
resummation at low px has been performed at the next-to-next-to-leading logarithmic level 

(nnll) nansinsi. 

Also at very high px, large logarithms arise, now of the form ln(l — px/Px 3X )i where p™ ax 
is the maximum kinematically possible transverse momentum for the vector boson at a given 
rapidity. For the photon p™ ax = , where Eqm is the machine center-of-mass energy (7 

TeV for the 2010-2011 LHC run) and y is the photon's rapidity. The approach to improving 
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on the fixed-order NLO calculation at high p T discussed in [T7J [TS1 [T9] was to expand around 
the limit p? = p™ ax . This is the machine threshold limit. When p? — £>™ ax , there is only phase 
space for the vector boson to be recoiling against a single parton, which also has pt = p™ ax . 
If the boson has slightly less pt, then the recoiling hadronic radiation must be jet-like, with 
the partons in the jet being either collinear or soft. Thus it is natural to describe the region 
near the machine threshold using Soft-Collinear Effective Theory (SCET) [20 | I2TI [22] . 

Using traditional methods, the threshold resummation for W/Z production at large pt 
was performed at NLL accuracy in [231 [2H [22]- Using SCET, the accuracy was increased 
to NNLL in [TSJ [19]. The effective theory approach simplifies the computations, and having 
operator definitions of the various ingredients of the factorization theorem lets us recycle 
known results, such as the 2-loop jet function, computed for other applications. This greatly 
reduces the amount of new analytical results needed. Nevertheless, also in the traditional 
formalism, the results were recently extended to NNLL accuracy [221 EZ] (although only the 
NNLO fixed-order expansion of the resummed result was computed in these papers). 

In practice, the threshold logarithms are important well away from the machine threshold 
because of the rapid fall off of the PDFs towards larger values of the momentum fraction x 
which ensures that most of the cross section comes from a region near the partonic thresh- 
old [2B1 [2H]- To what extent this dynamical enhancement of the threshold is effective was 
analyzed in detail in [30]. There is a simple phenomenological argument why it should hold: 
in events with a 300 GeV gauge boson, there is almost always a jet with pr ~ 300 GeV 
recoiling against it. That this jet is highly collimated and nearly massless indicates that the 
phase space region generating the large logarithms relevant for the vector boson px spectrum 
is important. Indeed, in cases such as inclusive Drell-Yan and Higgs production, where the 
NNLO corrections are known, it is found that 80 - 90% of the perturbative corrections to the 
cross section arise from the threshold terms, even in cases such as Higgs production, where 
the fall-off of the PDFs is not very strong. We thus expect our resummed results to provide 
a good approximation to the full NNLO result. 

2 Effective Field Theory approach 

The effective field theory allows us to obtain logarithmic contributions to the vector boson 
Pt spectrum which supplement the exact NLO distribution, computed in full QCD. These 
logarithmic terms arise from the threshold region, where the vector boson has the kinematically 
maximal transverse momentum. In this region, the jet recoiling against the vector boson is 
nearly massless. The formal derivation of these threshold terms is performed in the machine 
threshold limit, where x — > 1 for both PDFs. However, once the logarithms are extracted 
they can be used as additional information about the cross section in the kinematic region 
where x has more reasonable values. That is, the same threshold logarithms are present 
for any x, since they come from a perturbative calculation in QCD which factorizes from 
the non-perturbative PDFs. The only difference is that away from the machine threshold, 
we are no longer guaranteed that the threshold terms dominate the hadronic cross section 
parametrically. In practice they still give rise to the bulk of the cross section thanks to the 
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dynamical threshold enhancement discussed above. Physically, the threshold logarithms are 
associated with collinear radiation in the recoiling jet, or soft radiation coming from the jet or 
the incoming partons. A method to sum these logarithms to all orders in perturbation theory 
using the renormalization group (RG) in SCET was developed in J3TJ [32]. Its application to 
direct photon production was discussed in detail in [18]. Here we will only briefly summarize 
the method. 

A simple variable to use for the expansion near the machine threshold is Mx, the mass 
of everything-but-V^, where V refers to the vector boson (7, W or Z). In terms of the proton 
momenta and P% and the vector boson momentum q^, 

M 2 x = (P 1 + P 2 + q) 2 . (1) 

Since Pi and P2 are fixed, Mx is determined completely by the momentum of the vector 
boson. As its transverse momentum approaches its maximum allowed value at fixed rapidity, 
Mx — > 0. To understand the relevant degrees of freedom, it is helpful also to consider the 
partonic version of Mx, called mx- This is defined as 

m 2 x = (pi+P2 + q) 2 , (2) 

where p± = X\P^ and p% = X2P2 are the momenta of the partons coming out of the protons 
which participate in the hard interaction. Taking mx — > is called the partonic threshold 
limit. Obviously, Mx — > implies mx — > 0. The partonic mx is like Mx without including the 
beam remnants. Away from the machine threshold, the beam remnants make Mx large while 
mx can remain small. Thus the logarithms we actually expect to be important in affecting 
the vector boson px spectrum beyond NLO can be deduced by considering the theoretically 
simpler but less physical partonic threshold limit. 

Near the partonic threshold, the vector boson must be recoiling against a jet and there is 
only phase space for the jet to be nearly massless. So then we can write 

m 2 x = (pj + k s ) 2 ^p 2 j + 2Ejk, (3) 

where pj and kg are the collinear and soft momenta in the jet, Ej is the jet energy and 
k = pj ■ ks/Ej. It is because of this decomposition that the logarithmic terms we will extract 
come from either collinear effects (p 2 — > 0) or soft effects k — > 0. SCET implements the 
structure of the soft and collinear emissions on the Lagrangian level, using different fields 
to describe the soft and collinear partons. Via a field redefinition, the two sectors can be 
decoupled, after which the soft emissions are obtained from soft Wilson lines running along 
the directions of large momentum. 

The result from SCET is that the partonic cross section in the threshold region for any 
particular channel has the form 

s -, = a {0) (u,i)H(u,i,M v ,fj,) dk J(m 2 x - 2Ejk)S(k, y) , (4) 

dudt J 

where the partonic Mandelstam variables are s — {pi + P2) 2 , t — (pi — q) 2 and u — (p 2 — q) 2 , 
with q the vector boson momentum, with q 2 = My. We have factored out the Born level 
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Figure 1: Left: Factorization of the scattering amplitude near the partonic threshold. Right: 
Examples of NLO corrections to the hard, jet and soft function (from top to bottom). The 
thick blue lines denote partons collinear to the directions of the jet or the incoming hadrons. 
Soft emissions are pictured by thin red gluon lines. 

cross section a^(u,t). The hadronic cross section is obtained after convoluting with PDFs 
and summing over all partonic channels (see Sec. l2.1l below). 

The factorization theorem in Eq. (j4j) is depicted in Figure [U The hard function H contains 
the virtual corrections to the underlying hard-scattering process. There are two channels 
relevant for vector boson production, the Compton (qg — > Vq) and annihilation {qq — > Vg) 
channels, and the corresponding hard functions are related by crossing symmetry. A sample 
NLO contribution to the hard function in the annihilation channel is the top one-loop diagram 
on the right-hand side of Figure [TJ For the photon case, the one-loop hard function was given 
in [18], and in [19] it was outlined how the hard function can be obtained for my ^ 0. For 
completeness, we list the one-loop result for both the Compton and annihilation channel in 
the Appendix. The jet function J encodes the collinear emissions inside the final state jet, 
while collinear emissions along the initial state partons are absorbed into the PDFs. The jet 
function is obtained from the imaginary part of the two-point function of collinear fields (see 
the middle Feynman diagram on the right in Figured]). The two-loop results for the inclusive 
quark and gluon jet functions relevant here were obtained in [33] and [M]. The last Feynman 
diagram in the figure shows a NLO correction to the soft function, which describes the soft 
emissions from the energetic partons in both the initial and final state, which are encoded in 
Wilson lines along the corresponding directions. The corresponding soft function was recently 
computed to two loops in [35] . 

In the remainder of this section, we give the resummed result for the cross section and dis- 
cuss its numerical implementation. We first set up the integration over the parton momentum 
fractions in a form suited for threshold resummation and then give the resummed result, as 
well as the matching to fixed-order perturbation theory. Finally, we discuss how subtractions 
can be used to improve the convergence of the numerical integrations. 
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2.1 Integration variables 



The perturbative calculation, whether NLO or including resummation, produces partonic cross 
sections. The observable boson px spectrum is then obtained after convoluting with PDFs, 



dV 

dydpi 



ab 



^ / dZi / dx 2 f a /N 1 (x 1 , fi)f b/N2 {x 2 , fl) 



dydpx 



(5) 



where the sum is over all partonic channels, a, b G {g, g, (7}. The partonic cross section can 
also be written as 

d 2 cr ab A d 2 6- a6 



dydp 2 . didu 
At NLO, in a given channel, it has the general form 
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(7) 



where and the coefficients hi are functions of the two variables t and u. Because of the rela- 
tion s+t+u = m 2 x +My the 5-function parts, effectively, only depend on a single variable. The 
^-distributions are generalizations of the usual -(—distributions to dimensionful variables [36] . 
A N n LO computation would give distributions with logarithms up to ln 2n ~ 1 (mx//i) in the 
numerator. Resummation allows one to predict these singular terms at higher orders, but not 
the regular parts, such as . 

The leading-order cross sections for the production of a photon are 

2C F 7ra c . ra a s (/j 1 ) 
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Nr.8 



e T (u,t), 



a 



(0) 

qg 



7r a r , 



NJ 



-e q T {s,t) 



where e q is the charge of the quark and 
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(9) 



For the photon My 

M 



q 2 = 0, but we need the same expression also for Z and W bosons. 

(o) 



The amplitude ag q is obtained by replacing T (s,i) — > T (s,u) in the expression for a qg . To 
obtain the amplitude for Z production one replaces the quark charge in ([8]) by 

2 ls!l 2 + l^| 2 _ (l-2|e g |sin 2 ^) 2 + 4e 2 sin 4 ^ 



Ml 



8 sin 2 6 W cos 2 6 



(10) 



where 6w is the weak mixing angle. Since the W bosons have flavor- changing couplings, the 
sum over flavors must be replaced by a double sum over individual quark and antiquark flavors, 
q and q'. Only left-handed currents appear in this case. The relevant coupling for a W~ boson 
produced in the annihilation of an anti-up and a down quark is 



where V q > q are elements of the quark mixing matrix. 

Because of the singular behavior of the partonic cross section at threshold, it is advanta- 
geous to introduce m\ as an integration variable. Following [BJ, we perform the integrations 
in the form 

^^^.L^s + u-Mli m "d™i/W^)W^)|^, (12) 

with 

u = (P 2 - qf = M 2 V - yfiyjM*+p>eV , 

ml m = u + x 1 {M*-u), (13) 
—u 

^min = "T72 " 

M x — u 

When performing the resummation, one performs an expansion of the cross section around the 
partonic threshold mx = 0. With the choice of variables adopted in Eq. f fl2|) . the expansion 
is performed at fixed x± or, equivalently, at fixed t = (pi — q) 2 . This is problematic, since the 
expansion then induces unphysical rapidity asymmetries. In order to avoid this and obtain a 
symmetric form, we integrate twice: first with Eq. (fl2l) in the variables X\ and mx and then 
with the u <f-» t crossed version of Eq. (TT2|) in the variables x<i and nix- By taking the average 
of these two results, we obtain a symmetric form of the expansion around the threshold. In 
the case of direct photon production, a more convenient choice of integration variables is 

v = 1 + * , w = . (14) 

s s + t 

In this case, the partonic threshold is at w = 1. Using these variables significantly improves 
the numerical integration. The resummed photon cross section in v and w was given in [IS] . 

2.2 Resummation and matching to fixed order 

The resummed result for the cross section is obtained by solving the RG equations for the 
hard, jet and soft functions. Each function is then evaluated at its natural scale, where it 
does not suffer from large logarithmic corrections, and evolved to a common scale fi, which 
we identify with the factorization scale, see Figure |2j The solution of the RG equation for the 
hard function in the annihilation channel is 

cr Q qq [u, t, fi) H qq (u, i, fi) = U^(ii h , fi) a° qi! (u, i, ii h )H qq (u, i, fi h ) . (15) 

The matching scale \Xh is chosen to be of the order of pr to avoid large logarithms. The 
NLO result for H qq (u,t, fih) was given in [19]. We have included the Born cross section which 
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Figure 2: Resummation by RG evolution. 



depends on the scale via the coupling constant a s (fi). The evolution factor 11^(^,1^) for the 
above combination takes the form 



\nU^(jjL h ,tj) = 2(c F + 



2S{n h ,n) - A cnsp (n hl fi) In -j 
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-2^(^,//), (16) 



"^cZa^^. (17) 
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Explicit expressions for these functions in RG-improved perturbation theory can be found 
in [32]. The function is the same as A cwsp (v, //) with 

Ca s 2 

1h cM = 27, + 7 3 - — In 7cusp (18) 

replacing 7cU sp- The quark and gluon anomalous dimensions 7 g and 7 g are given to three-loop 
order in [3Tj . The evolution factor [/^ can be obtained from the above results using the 
crossing relation s -h- — u at fixed t , and t/ji follows from using £ — > u. The resummed 
results for the jet and soft functions can be obtained by solving their RG equations in Laplace 
space [31]. For the gluon jet function, for example, the result takes the form 



1 ip 

J 9 (p 2 ,/i) = U jg (fx j , f i)j g {d Vjg )— -2 



2 \ %s e -7EVjg 



where j g is the Laplace transform of the momentum-space jet function and 

UjgiP-vl*) = exp[-4C A S'(// i ,//) + 2A Jg (p j ,fj,)] , 
Vj g = 2C Y A^4cusp(A t j; aO • 



(19) 
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The corresponding results for the quark jet function and the soft functions in the different 
channels are all listed in [T8], together with the necessary anomalous dimensions. Inserting 
the resummed expressions into P]h one obtains 

d 2 <r - 

-2L = a®jj fi h ) H q g(u,i, fi h ) U Hq9 (n h , //) U Jg (nj, fi)U Sq9 (n s , fi) 



' J m 2 x \ fij J \E h fi s J r^j + rjs) 

a${u, i, n h ) H q q(u, i, fi h ) U Hq9 (y h , fi) U Jg (^, y) U ScM {fi s , y) 

x i 9 K ? , 3^ + In j —3 , (20) 



where E h = \Jiu/As = pr/2 and 

Vqq = r fg + Vgg = 2CUA cusp (/ij, /i) + (4C F - 2C A )A cnsp (fi s , (J>) ■ 

To arrive at this expression, the convolution integral in fl4]) was explicitly carried out, which 
is possible because of the simple form of the RG evolved soft and jet functions. 

Using the general expression Eq. ( 120]) for the resummed cross section, we can now explicitly 
compute the resummed distribution. We include almost all of the ingredients for N 3 LL accu- 
racy. For N 3 LL, one needs the cusp anomalous dimension to four-loop order and the regular 
anomalous dimensions to three loops, together with the two-loop results for the hard, jet and 
soft functions. For the anomalous dimensions, the only missing ingredient is the unknown 
four-loop cusp anomalous dimension, which we estimate using the standard Pade approxima- 
tion as T4 = r3/(r2) 2 . The second ingredient, which we do not include are the non-logarithmic 
pieces of the jet, soft, and hard functions. The full two-loop jet functions are known [33j El] 
and also the two- loop soft function has now been computed [35] . The non-logarithmic piece of 
the two- loop hard function can be extracted from the results of J38J [39] , and we plan to include 
the full two-loop matching in the future. To indicate that we only have a partial result, we 
denote our highest order by N 3 LL p . 

In order to include all the available perturbative information, we match to the NLO fixed- 
order result, which is the highest order available. To perform the matching, we use 



/ ,2 \ N 3 LL p +NLO / jo x N 3 LL P / , 2 \ NLO / , 2 \ NNLL 

/ d z a \ l d z a \ ( d z a \ ( d l ~ 



i era \ / a~a \ r / a~cr \ / a-<r 

Kdydp^J \dydp^J \dydp^J \dydpj 



(21) 



The subscript on the last term indicates that all scales must be set equal to the relevant value 
of /i. Setting these scales equal switches off the resummation. The NNLL expression includes 
the one-loop corrections to the hard, jet and soft functions. Once it is evaluated with all scales 
equal, it reduces to the singular threshold terms of the NLO result, which must be subtracted 
since they are already included in the resummed result. 
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2.3 Subtractions 



To compare with data, one needs to perform 4-dimensional integrals, over xi, mx, V and pr- 
These numerical integrals are computationally expensive and additionally challenging because 
of the singular nature of the partonic cross sections. After resummation the partonic cross 
sections are no longer distribution valued at the partonic threshold, but behave as 



m\ \ /i 2 



m 



x 



(22) 



near the threshold, see (1201) . For the natural hierarchy fij > /j, s > [i the quantity r\ is larger 
than zero and the integral over mx converges. To see this one rewrites r) qq - in the form 

r] q g = 2C A A cusp (fij, /i s ) + AC F A cusp (fi s , ji) . (23) 

In practice, however, the scale hierarchy is not very large, so that convergence can be quite 
slow. Furthermore, we will choose a high value of the factorization scale /i, in which case the 
integral is no longer guaranteed to exist since r] qq - can become negative. 

For some threshold variable m, the integral we need to evaluate has the form 

J(M) = / dmm^/H. (24) 



o 



To analytically continue the result to negative r\ values, or to improve convergence, it is useful 
to perform subtractions [301 E2] • A single subtraction would use 



J(M) = J" dmm^j/ffl) + [/(m) - /(0)] } (25) 

M r i i 

dm j-M 2 "- 1 /^) + m 2 "- 1 [/(m) - f(0)} j , (26) 

where the difference f(m) — f(0) makes the integral more convergent. Indeed, assuming f(m) 
is smooth the integral will now converge for rj > —1/2. A second subtraction would give 

J(M) = ^ M dm|^M 2 ^V(0) + ^^M 2 V'(0)+m^-M/(m)-/(0)-mr(0)]| , (27) 

which makes the integral converge for rj > — 1. Analogously, one can perform higher-order 
subtractions to make the integral more and more convergent. In practice, performing too 
many subtractions slows the code down because the expressions become lengthy and the 
subtraction itself can become numerically unstable. Generally, we have found that using two 
or three subtractions gives stable results for the N 3 LL resummed integrand. Note that higher 
subtractions involve derivatives of PDFs, which we compute by first interpolating the PDFs. 
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3 Scale setting 



The resummed distribution obtained using SCET involves four matching scales: the hard scale 
Hh, the jet scale fij, the soft scale /i s and the the factorization scale fi. Each ingredient of 
the factorization formula, thus, can be evaluated at its appropriate scale. This is in contrast 
to the fixed-order calculation which involves only a single scale, the factorization scale /i. In 
addition to the factorization scale, a renormalization scale fi r is often introduced by hand in 
fixed-order computations, by expressing the coupling constant as a s (fi r ). Introducing a second 
scale n r may be useful as a tool to estimate uncertainties, but there is no physical justification 
for having \i ^ fi r when working at fixed order. In contrast, the additional scales in SCET 
correspond to different physical regions. 

If one chooses all the scales equal in the SCET calculation, the resummation is switched off. 
Then only a single renormalization group scale remains. In this limit, SCET generates all of 
the terms that are singular in the threshold limit at next-to-leading order, as shown in Eq. ([7]). 
Since we include all logarithmic pieces in the two-loop hard, jet and soft functions, we also 
obtain all singular terms at NNLO, except for the coefficient of 8{m 2 x ). In general, the singular 
terms amount to a large fraction of the full perturbative correction. For electroweak boson 
production at large transverse momentum, they amount to 70-80% of the NLO correction |19j . 
For inclusive Drell-Yan and Higgs production, also the NNLO correction is known, and it is 
found that a similarly large fraction of the perturbative corrections arises from the partonic 
threshold region [301 SO] • 

A perpetual frustration with fixed-order calculations is that they provide no insight into 
how to choose the factorization scale in problems that involve physics at several scales. For 
example, for W production either jj, = mw or fi — \J + p\ might seem natural. For large 
Pti the difference in the prediction between these different parameterizations is larger than 
the scale variation within any particular parametrization. See for example Figure El to be 
discussed more below. If higher order calculations were available, as they might be soon for 
direct photon or W production, the scale variation and parametrization dependence would 
weaken. Unfortunately, however, there are only a handful of observables which have ever 
been computed beyond NLO. For more complicated processes, such as W+4 jet production, 
NNLO is a distant hope, and in this case there are many possible natural parameterizations. 
Thus the choice of parametrization for the factorization scale can amount to a significant and 
difficult-to-estimate source of uncertainty for a fixed-order computation. 

An extremely satisfying feature of the effective field theory approach is that it does indicate 
what the appropriate parametrization should be. For some observables, such as e + e~ event 
shapes [JTJ H2J H3j 01] the scales are manifest in the resummed distribution. In hadronic 
collisions, one needs a somewhat more sophisticated procedure since the scales can depend 
on the functional form of the non-perturbative PDFs. A method for determining these scales 
without any arbitrary input of what is natural and what is not was proposed in [30] and applied 
to direct photon and W production in [181 US]- The idea is very simple: one supplements the 
leading-order calculation with just one part of the SCET calculation at a time, for example, 
the hard, jet or soft function evaluated at NLO in fixed-order perturbation theory. Doing this, 
there should be a single scale (p) which is the average value of momentum associated with 
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Figure 3: Scale sensitivities. These plots show the effect of adding part of the fixed-order NLO 
calculation to the LO calculation. The left panel shows what happens if all the /i-dependent 
terms at NLO are added together. There is a slow monotonic logarithmic [i dependence, 
with no natural extremum. In contrast, when the hard, jet, or soft contributions are added 
separately, there are natural extrema. These extrema indicate the average value of momenta 
(p) appearing in the logarithms. That there are different extrema for the different components 
proves that multiple scales are relevant. The plots are for W + bosons, but the qualitative 
features are the same for all bosons. 



these degrees of freedom, appearing in the large logarithms. After integrating over the PDFs, 
the perturbative correction will then have the form 

4j = a s (fi)(c 2 L 2 + Cl L + c ) , (28) 

with L = l n 7^y- If /^ is chosen either much lower or much higher than (p), the perturbative 
corrections will become large. Since we do not have an analytic expression for the distribution, 
due to the necessity of convoluting with PDFs, we determine (p) numerically by computing the 
individual corrections to the cross section as a function of [i. The result is shown in the right 
plot of Figure [3J It has the expected form ([25]) and we see that while the jet and soft scales 
are concave upwards, the hard curve is concave downward. The extrema of the corresponding 
curves indicate the scales (p) that dominate these contributions after integrating over the 
PDF. It is then natural to define our default values for /j, as the positions of the extrema. 
That there are different extrema for the different components proves that multiple scales are 
relevant. These scales are conflated in the fixed-order calculation. The left plot in Figure |3] 
shows the fixed-order scale dependence. In this case, there is monotonic [i dependence, with 
no natural extremum. 

To find the scales numerically we extract these extrema from the curves. Using a number of 
different machine center-of-mass energies (we tried 2, 7, 14, and 100 TeV), pp and pp collisions, 
and various boson masses, we determine a reasonable approximation to these points is given 
by the following functional forms 
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Figure 4: Natural hard and jet scales. These scales are determined, after integrating over the 
PDFs, as typical energies appearing in the logarithms when the hard or jet functions at next- 
to-leading order only are included. Dots show extrema for various values of boson masses, 
and lines show our approximations. An important qualitative point is that the jet scale is 
naturally lower than the hard scale. This is an output from our numerical procedure, not an 
input from a formal analysis. 



13p T + 2M V p 2 T 

» h = 12 -TS' (9) 

7Pt + 2M V f 2p T \ 

N = — To — 1 - ^ • ( 30 ) 



12 V 

A comparison of this fit to the extrema for the hard and jet scales is shown in Figure HI We 
have constrained the jet scale to vanish at the endpoint pt = y/s/2 since at that point there 
is no phase space for collinear emission and the recoiling jet must be massless. An alternative 
and slightly simpler hard scale choice that is quantitatively equivalent for LHC energies is 

7PT + My 

Vh= g • (31) 

In the comparison to data, we use the scales in Eq. ff29l) . for consistency with [19] . which used 
these scales in comparison to Tevatron data from run II. 

From Figures |3] or [H it is obvious by eye that the natural jet scale is lower than the 
natural hard scale. While the hard scale is actually fairly close to a common scale choice 
\i = \/m v + pj,, as shown on the left of Figure [51 the jet scale is significantly lower. Both 
scales are higher than the fixed scale /x = my which was used in comparison to mcfm in the 
atlas study of the W spectrum [2]. 

To further emphasize the importance of scale choices, we show in Figure the relative 
difference in the NLO prediction from the different parameterizations for the W spectrum. 
The band corresponds to a region |/x(pt) < A* < 2[i(pt) where fi(pr) is either mw, our hard 
scale, or the popular choice \ffn%r + Pt- For pt ~ m w, all scales give comparable results while 
for large pt, the fixed scale gives a prediction significantly higher than either of the other two 
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Pr(GeV) Pr(GeV) 

Figure 5: Comparison of our scale choices (blue) with the traditional choice /i = a/ m 2 ^ + p\ 
(red) for photon (dashed) and W boson (solid). For the hard scale, there is not much difference. 
On the other hand, the natural jet (and soft) scales are lower than the traditional choice, but 
higher than the fixed scale \i = my/. 

parameterizations. If we used the jet scale instead of the hard scale, the band would be closer 
to the fi = my/ band. Thus it is important to choose the appropriate scale in the appropriate 
place, to get an accurate prediction. 

Having determined the default values for the scales in Eq. (1291) . we can compute the 
resummed distribution. As discussed above, we include all ingredients for N 3 LL accuracy, 
except for the two-loop non-logarithmic terms in the hard, jet and soft functions. We match 
to NLO fixed order and denote our highest order resummed result by N 3 LL p +NLO, where 
the subscript "p" stands for partial. Convergence in perturbation theory and the relative 
size of various scale variations are shown in Figured To generate the bands in this plot, we 
determine the maximum and minimum cross section obtained when varying each scale up and 
down by a factor of 2 around its default value. In contrast to the fixed-order result, the scale 
dependence is not monotonic (cf. Figure [3]). To determine the maximum and minimum, we 
compute the cross section at |, 2 and the central value, fit a parabola to those three points 
and take the maximum and minimum along the parabola. 

Curves in the first four panels of Figure [7] are not matched to fixed order. The relatively 
large factorization scale uncertainty comes about because the /i dependence is only canceled 
in the resummed distribution near threshold. The full \x dependence at NLO is removed once 
the theory is matched to the fixed-order distribution, as it is in the bottom two panels. The 
combined uncertainty that we use for our final error estimates is the quadratic sum of the 
hard, jet, soft and factorization scale uncertainties: 

Aa = yj (A h a) 2 + (A j( x) 2 + (A s cr) 2 + (A /C r) 2 (32) 

This is a conservative error estimate. We observe that the N 2 LL and N 3 LL P scale variation 
bands overlap, but the increase in the cross section from NLL to N 2 LL is larger than the 
NLL band. The same behavior is also seen when going from the LO to the NLO fixed-order 
result. The corresponding bands would overlap had we evaluated the LO and NLL results with 
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Figure 6: Scale variations at next-to-leading order. The blue southeast stripes show the scale 
variation of the NLO calculation (called NNLO in fewz) with n = [if = \i r = mw, as in the 
atlas paper. The red northeast stripes show the prediction using /j,f — fx r — \J m w + p^, and 
the black vertical stripes have fif and \x r set to the scales in Eq. ( 129]) . Bands correspond to 
varying /i = /ij = /i r by factors of two from these default scales. 

leading-order PDF sets which have a larger value of a s , instead of the NNLO PDFs we use 
throughout. The increase in the cross section from NLL to N 2 LL is mostly due to the one-loop 
constants in the soft and hard functions, as can be seen from the right panel of Figure [3j We 
have checked how much of a shift the known two-loop jet and soft function constants induce 
and find that it is below a per cent. 

4 Comparison with LHC data 

We are now ready to compare to LHC data. We discuss separately the two processes we study, 
direct photon and W production. For numerical work we use the NNLO MSTW 2008 PDF set 
and its associated a s {m z ) = 0.1171 jl5]. We also use mw = 80.399 GeV, a c . m . = 127.916" 1 , 
sin 2 ^ = 0.2226, V ud = 0.97425, V us = 0.22543, V ub = 0.00354, V cd = 0.22529, V cs = 0.97342 
and = 0.04128. 

4.1 Direct photon 

For direct photon production, to be consistent with the comparison to Tevatron data in [18] , 
we use the scale choices from that paper 

fJ-h = Pt , 
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Figure 7: The cross section for W + production at the LHC 7 TeV for various orders in 
perturbation theory, normalized relative to the next-to-leading fixed-order curves. In the top 
four panels the resummed curves are not matched to fixed order, which shows how including 
just the logarithms compares to the full result. Bands come from variations of the hard, 
jet, soft and factorization scales by factors of 2 around our default scales are shown and 
taking the maxima within that variation region. The fifth panel shows the factorization scale 
variation after matching. The sixth panel is the uncertainty from adding the hard, jet, soft 
and factorization uncertainties in quadrature. 
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These are slightly simpler than those in Eq. f[2"9~j) . but equivalent within the uncertainties. 

The direct photon spectrum is complicated by the requirement of photon isolation which is 
necessary to remove hadronic backgrounds, such as n° decays. The atlas study [TJ required 
energy in an cone of R = 0.4 around the photon to have energy less than 4 GeV. To include 
this effect rather than matching to the inclusive fixed-order NLO calculation, we match to 
the NLO calculation with isolation and fragmentation contributions using the Monte Carlo 
program JETPHOX [4"6] . 

The atlas data in [1] includes 35 pb" 1 of data separated into three rapidity regions. The 
comparison of the theory prediction at NLO and at N 3 LL p +NLO order is shown in Figure [HJ 
The theory and data are in agreement within uncertainties. It will be interesting to update 
this comparison to a large data set, particularly if it includes higher p T . 

4.2 W boson 

The calculation of the W boson px spectrum at N 3 LL p +NLO order is significantly more 
challenging numerically than direct photon production, despite the identical factorization for- 
mulae. The extra scale, the boson mass, complicates the kinematics, which makes the integrals 
converge more slowly and the scale choices more complicated. Moreover, experimentally, since 
only the charged lepton from the W decay is measured, acceptances have to be included in 
comparing the inclusive W spectrum to the measured distribution. 
The acceptance cuts used by atlas in [2] were 

p l T > 20 GeV, \rf\ < 2.4, p v T > 25 GeV , (34) 

and 

= \/ {\P\\ + \p^\) 2 - \pV + Pt u \ 2 > 40 GeV. (35) 

where / is an electron or muon. 

To apply these cuts to our inclusive sample, we multiply our inclusive cross sections by 
acceptances, given by the ratio of the inclusive W cross section to the W cross section in 
the fiducial volume as a function of px- These acceptances are shown in Figure EO We 
calculate them at leading order and next-to-leading order for W + and W~ using the program 
fewz [12]. We find no significant difference between the leading order and next-to- leading 
order acceptances. For the numerical work, we use a smooth function fit to the inclusive NLO 
acceptances for W + + W~ , also shown in Figure M 

The total cross sections we find at the -y/i = 7 TeV LHC, using jjf = [i r = mw are 



a(pp -> W + -> fi + u) inc = (6204 ± 0.7) pb , 
a{pp W~ -»■ ^u) inc = (4326 ± 0.5) pb , 



a{pp ->■ W + -> fi + u) M = (3061 ± 3.0) pb , 

a{pp -»■ W~ -»■ /i~z/) fid = (2038 ± 1.9) pb , 

(36) 
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Figure 8: Comparison of theory to ATLAS data for direct photon in various rapidity regions. 
The normalization is such that corresponds to the inclusive direct photon spectrum at NLO 
without isolation or fragmentation. The NLO cross sections including fragmentation and 
isolation are slightly higher, especially at low pt, since the enhancement due to fragmentation 
outweighs the reduction due to isolation. Uncertainties on the NLO correspond to factor 
of 2 variation of /!/ = fi r around its default value of pr- For the resummed curves, which 
also include isolation and fragmentation, uncertainties are hard, jet, soft and factorization 
variations by factors of 2 added in quadrature. Uncertainties from fragmentation and isolation 
are not included, and may be large. 
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Figure 9: Kinematic acceptances according to atlas cuts. There is little difference in the 
acceptances for W + or W~ bosons. However, there is about a 2% increase in the efficiency 
when going from leading-order to next-to-leading order in perturbation theory (all computed 
with fewz). The NLO acceptances are used in our comparison with atlas data. 

where "inc" refers to the inclusive cross section and "fid" to the fiducial cross section (with 
cuts). Errors are integration errors from FEWZ. Dividing these by the branching ratio to 
muons, BR{W — > p,v) = 0.1083 gives the total inclusive cross section. To compare to data, 
we take our theoretical calculation of the inclusive differential cross section, multiply by the 
acceptance curves and then divide by the total cross section in the fiducial volume, 47.04 nb. 
This lets us compare directly to the atlas data, which is normalized to the total number of 
events in the fiducial volume. 

The comparison to the ATLAS data is shown in Figure [TDJ The agreement is excellent. In 
this plot, results are shown normalized to the NLO prediction with \i$ — /j, r = raw- As we 
have argued in the previous section this is not a good scale choice in the large pr region. We 
use \i = mw as the basis for our comparison since it is the scale choice used by atlas in [2], 
and, therefore, our NLO calculation can be directly compared to their calculation labeled 
MCFM in Figure 7 of [2]. 

There are two important qualitative conclusions that can be drawn from Figure [10J The 
first is that scale choices are important. Although the data is within the NLO uncertainty, 
comparing the NLO band in this plot, which uses a fixed scale \x = my/ with those in Figure El 
we see that the downward trend in the data indicates a clear preference that /x should increase 
with pt- However, as we have argued, there is no natural scale choice at fixed order, since 
multiple scales are present. Our procedure of determining these scales numerically is supported 
quantitatively by the agreement between the N 3 LL p +NLO band and the data in this figure. 

The second qualitative conclusion concerns the PDF uncertainties. These are shown as the 
blue dashed lines in Figure [TD1 The PDF uncertainties are smaller than the uncertainties on 
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Figure 10: Comparison of theory to atlas data for the W spectra. The red band is the NLO 
prediction, using fif = \i T = m w , as in [2]. The N 3 LL p + NLO prediction, in green, is in 
excellent agreement with the data. Dashed blue lines indicate PDF uncertainties which are of 
order the scale uncertainties at N 3 LL P + NLO order. 

the data and the NLO scale uncertainties, but of the same order as the scale uncertainties of the 
resummed distribution. This indicates that PDF fits could be improved using the W spectra, 
but only if resummation is included (or perhaps if the NNLO result becomes available). 



5 Conclusions 

In this paper, we have compared theoretical predictions for the direct photon and W boson 
spectra at high p? to measurements performed by the ATLAS collaboration using LHC data. 
The predictions were performed using the exact cross section at next-to-leading order in a s 
(the highest order known), supplemented with additional terms to all orders in a s coming 
from a threshold expansion. These extra terms correspond to large logarithms associated with 
infrared singularities of the recoiling jet. To isolate these terms, the resummed calculation is 
performed near the partonic threshold, in which the Pt of the vector boson is maximal for the 
given value of the partonic center-of-mass energy. In this limit, the cross section factorizes 
and the logarithmic terms can then included to all orders in perturbation theory. These 
terms usually give the dominant contribution to the the cross section. In the photon case, 
the fragmentation cross section and isolation corrections were added also, using the program 

JETPHOX. 

A main advantage of the resummed cross section, which was calculated using effective 
theory in [18] and [19], is that it has well-defined scales associated with different phase space 
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regions. Unlike a fixed-order calculation, which merges all the scales into one, the effective 
field theory allows one to choose the scales appropriate for the relevant regions. We employed 
a numerical procedure to determine which scales are appropriate. This removes a source of 
uncertainty from the fixed-order calculation, namely, which value of the scale /i should be 
adopted. Typical fixed-order calculations choose scales like mw, Pt or a/Pt + m w- We find 
that a single such parametrization is insufficient: while the hard scale is naturally close to 
\/p\ + mw, the jet and soft scales are naturally lower. 

The results of our comparison with data are shown in Figures IBl and fTUl The photon case 
is complicated by the requirement of photon isolation, while the W case is complicated by the 
missing energy and the necessary acceptance cuts on the lepton. We found good agreement for 
the direct photon case, and excellent agreement for the case of the W boson. In the W boson 
case, one can clearly see the importance of proper scale choices. Moreover, the reduction of 
theoretical uncertainty when resummation is included is enough to make it comparable to the 
PDF uncertainty. These comparisons provide a convincing demonstration of the relevance of 
resummation for LHC physics. 

In the future, it would be useful to compare our prediction to the direct photon data 
available from cms [UJ . It would also be interesting to compare to Drell-Yan spectra from 
intermediate Z bosons at high pt[1HI HH], and to direct photon and W data at higher lumi- 
nosity. Improvements of the theoretical description could be achieved by including the full 
two-loop hard, jet and soft functions and the fixed-order NNLO calculation, once it becomes 
available. Furthermore, on top of the QCD effects, one should also include electroweak Su- 
dakov logarithms [50~ll5T]. which will have a noticeable effect at high p?. In addition, we hope 
to eventually provide a publicly available code to produce the resummed results. 
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A One-loop hard function 

In this appendix we give the result for the hard function in both the annihilation and the 
Compton channel. The functions are related by crossing symmetry, but the analytic continu- 
ation from one channel to the other is not entirely trivial. The hard function is obtained from 
the result (A. 9) in [7J after performing renormalization. 
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For the annihilation channel, we have 



H g q(u,t) = 1 + 



Air 



7T 



6 



-16 



7tt 2 



2C A ln 2 -^ + C A ln 2 ±v±v 



Ml - t 



Ml 



9 Ml — u , s 
+C A In 2 y „ + In 



M 2 



M 2 



s 2 \ 



-6C F - 2C A In — - C A In 



M 2 



- 6C F In 



-2C A In — In 
tu 



-C A - 2C F ) In 



2 



+2CU Li s 



My 



Ml-t 



2C A Li2 



My 



My - U 



+ 



{-C A + 2CV) 
t 



c f 



s + t 



+ [C A 

+ (c A 



s + u 
u 



+ c, 



s + t u 

Ml (t 2 + u 2 ) 
tu (t + u) 

4s 2 + 2st + 4su + tu 



s s + u 
+ — 



s + u 



(t + u y 



2s 
t + u 



In 



Ml 



In 



M 2 



(s + uY 
As 2 + 4st + 2su + tu\ , -u 



s + t 
C A + 2C F ) 



(s + t) 

2 



hi 



Ml 



s + (s + u 
tu 



-ln 2 ^, 

2 My 2 



1 . 2 M 2 - t 



-In 



M 2 



+ In 



In 



-t 



Ml s-M 2 



Lin 



My 



M 2 V 
Ml-t 



s 2 + (s + ty 

tu 



+ Li 2 



M 2 V 



-hv 

2 

Li 2 



-u 



_^_1 2 Ml-u ln _f_ 
M 2 2 M 2 f M s-M 



M 2 



My — M 



and the Compton channel result reads 



a 



H qg (u,t) = l + -^-lC A — + C F 



4:71 



7tt 2 



6 



-16 



7T 



„ , 2 — St 



6C F In — ft — CU In 

M 2 Mi 



+ C A hv 



M 



-2Cln^i^ln^ + C,ln 2 



MyW 



M 2 



M 2 



" 2^1n (M ^7t )gt ln- M 



MyM 2 



-2C A In 2 2C F In 2 — + -6CV + 2C A In - + 4C F In 

u 



-u 



In 



21 



-(C A + 2C F ) \n 2 ^-2C A U 2 

s \ s 



M 2 \ / M 2 



M 2 -t 



+ 



4vr ^(°) 
C A 



s 



C f 



it 



S + U 



s + u 

t 



u 



t + u 



t + u 



t + u 

t 



' s + u 
+ (-C A + 2C F ) 

- (-C A + 2C F ) 



st + 2su + Atu + Au 2 
(t + uf 

st + Asu + 2tu + 4m 2 
? (s + uf 

M 2 (s 2 + t 2 ) | 
st (s + t) 



111 



111 



M 2 



-t 



2u 



u 



s + t (s + ty 



In 



-u 



M 2 



U 2 + (t + Uf (1 ^ S 



+ Li 2 



M 



Lie 



st 



M 2 
My — n 



2 My 2 



M 2 V 



Ml s-Ml 



U 2 +(S+U) f IT 2 1 o My — t 1 n My — M , -i, — It 

H — In 2 — 5-= In 2 — ¥—= h In — T In — T 

st V 2 2 M 2 2 M 2 M 2 M 2 



Li 2 



My 



Ml-t 



14' Mf 



M 2 -u 



(3* 



References 



[1] G. Aad et a/. [ATLAS Collaboration], Phys. Lett. B 706, 150 (2011) [ arXiv:1108.0"253l 
[hep-ex]]. 



[2] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 85, 012005 (2012) [arXiv: 1108.63081 
[hep-ex]]. 

[3] P. Aurenche, A. Douiri, R. Baier, M. Fontannaz and D. Schiff, Phys. Lett. B 140, 87 
(1984). 

[4] P. Aurenche, R. Baier, M. Fontannaz and D. Schiff, Nucl. Phys. B 297, 661 (1988). 

[5] L. E. Gordon and W. Vogelsang, Phys. Rev. D 48, 3136 (1993). 

[6] R. K. Ellis, G. Martinelli, R. Petronzio, Nucl. Phys. B211, 106 (1983). 

[7] P. B. Arnold, M. H. Reno, Nucl. Phys. B319, 37 (1989). 

[8] R. J. Gonsalves, J. Pawlowski, C. -F. Wai, Phys. Rev. D40, 2245 (1989). 



22 



[9] R. Gonsalves, http://www.physics.buffalo.edu/gonsalves/ 
[10 



[11 
[12 

[13 

[14 
[15 

[16 

[17; 

[18 
[19 

[20 
[21 
[22 
[23 
[24 
[25 
[26 
[27 
[28 
[29 
[30 
[31 



J. Campbell, K. Ellis, C. Williams, http://mcfm.fnal.gov/ 



K. Melnikov and F. Petriello, Phys. Rev. D 74, 114017 (2006) [hep-ph/0609070] . 

R. Gavin, Y. Li, F. Petriello and S. Quackenbush, Comput. Phys. Commun. 182, 2388 
(2011) [arXiv:1011.3540l [hep-ph]]. 

S. Catani, L. Cieri, G. Ferrera, D. de Florian and M. Grazzini, Phys. Rev. Lett. 103, 
082001 (2009) |arXiv:0903.2T20l [hep-ph]]. 

C. Balazs and C. P. Yuan, Phys. Rev. D 56, 5558 (1997) [hep-ph/9704258| . 

G. Bozzi, S. Catani, G. Ferrera, D. de Florian and M. Grazzini, Phys. Lett. B 696, 207 
(2011) [arXiv: 1007.23511 [hep-ph]]. 

T. Becher and M. Neubert, Eur. Phys. J. C 71, 1665 (2011) [arXiv:1007.4"005l [hep-ph]]. 

E. Laenen, G . Oderda and G. F. Sterman, Phys. Lett. B 438, 173 (1998) 
|hep-ph/9806"467| . 

T. Becher, M. D. Schwartz, JHEP 1002, 040 (2010). 

T. Becher, C. Lorentzen and M. D. Schwartz, Phys. Rev. Lett. 108, 012001 (2012) 
[arXiv:1106.4"3l0l [hep-ph]]. 

C. W. Bauer, S. Fleming, D. Pirjol and I. W. Stewart, Phys. Rev. D 63, 114020 (2001). 

C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. D 65, 054022 (2002) 

M. Beneke, A. P. Chapovsky, M. Diehl and T. Feldmann, Nucl. Phys. B 643, 431 (2002) 
N. Kidonakis, V. Del Duca, Phys. Lett. B480, 87-96 (2000). 
N. Kidonakis, A. Sabio Vera, JHEP 0402, 027 (2004). 

R. J. Gonsalves, N. Kidonakis, A. Sabio Vera, Phys. Rev. Lett. 95, 222001 (2005). 
N. Kidonakis and R. J. Gonsalves. larXiv: 1201 .52651 [hep-ph]. 
N. Kidonakis and R. J. Gonsalves. larXiv:1109.281~7l [hep-ph]. 

D. Appell, G. F. Sterman and P. B. Mackenzie, Nucl. Phys. B 309, 259 (1988). 



S. Catani, M. L. Mangano and P. Nason, JHEP 9807, 024 (1998) [hep-ph/9806484| . 
T. Becher, M. Neubert, G. Xu, JHEP 0807, 030 (2008). 

T. Becher, M. Neubert, Phys. Rev. Lett. 97, 082001 (2006). |hep-ph/0605'050] . 



23 



[32] T. Becher, M. Neubert and B. D. Pecjak, JHEP 0701, 076 (2007) |hep-ph/0607228] . 
[33] T. Becher, M. Neubert, Phys. Lett. B637, 251-259 (2006). |hep-ph/0603l"40] . 
[34] T. Becher, G. Bell, Phys. Lett. B695, 252-258 (2011). 

[35] T. Becher, G. Bell and S. Marti, JHEP 1204, 034 (2012) |arXiv:1201.5"572l [hep-ph]]. 

[36] F. De Fazio and M. Neubert, JHEP 9906, 017 (1999) |arXiv:hep-ph/9905351] . 

[37] T. Becher and M. Neubert, JHEP 0906, 081 (2009) [arXiv:0903. 11261 [hep-ph]]. 

[38] L. W. Garland, T. Gehrmann, E. W. N. Glover, A. Koukoutsakis and E. Remiddi, Nucl. 
Phys. B 642, 227 (2002) |hep-ph/0206067] . 

[39] T. Gehrmann and L. Tancredi, JHEP 1202, 004 (2012) jarXiv: 1112.15311 [hep-ph]]. 

[40] V. Ahrens, T. Becher, M. Neubert and L. L. Yang, Eur. Phys. J. C 62, 333 (2009) 
[arXiv:0809.4283l [hep-ph]]. 

[41] M. D. Schwartz, Phys. Rev. D 77, 014026 (2008) |arXiv:0709.2"709l [hep-ph]]. 

[42] T. Becher and M. D. Schwartz, JHEP 0807, 034 (2008) |arXiv:0803.0"342l [hep-ph]]. 

[43] Y. -T. Chien and M. D. Schwartz, JHEP 1008, 058 (2010) [arXiv: 1005. 16441 [hep-ph]]. 

[44] I. Feige, M. Schwartz, I. Stewart and J. Thaler. larXiv: 1204.38981 [hep-ph]. 

[45] A. D. Martin, W. J. Stirling, R. S. Thorne, G. Watt, Eur. Phys. J. C63, 189-285 (2009). 

[46] S. Catani, M. Fontannaz, J. P. Guillet and E. Pilon, JHEP 0205, 028 (2002) 
[hep-ph/0204023| . 

[47] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. D 84, 052011 (2011) 
[arXiv: 1108.2041 [hep-ex]]. 

[48] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. D 85, 032002 (2012) 
[arXiv:1110.4973l [hep-ex]]. 

[49] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 705, 415 (2011) [ arXiv:1107.2"38l1 
[hep-ex]]. 

[50] J. H. Kiihn, A. Kulesza, S. Pozzorini, M. Schulze, Phys. Lett. B609, 277-285 (2005); 
Nucl. Phys. B727, 368-394 (2005); Phys. Lett. B651, 160-165 (2007); Nucl. Phys. B797, 
27-77 (2008). 

[51] W. Hollik, T. Kasprzik, B. A. Kniehl, Nucl. Phys. B790 (2008) 138-159. 



24 



